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Abstract 

We extend our earlier work, regarding the perturbative stability of string configurations 
used for computing the interaction potential of heavy quarks within the gauge/gravity 
correspondence, to cover a more general class of gravity duals. We provide results, 
mostly based on analytic methods and corroborated by numerical calculations, which 
apply to strings in a general class of backgrounds that encompass boosted, spinning 
and marginally-deformed D3-brane backgrounds. For the case of spinning branes we 
demonstrate in a few examples that perturbative stability of strings may require strong 
conditions complementing those following by thermodynamic stability of the dual field 
theories. For marginally-deformed backgrounds, we find that even in the conformal case 
stability requires an upper value for the imaginary part o of the deformation parameter, 
whereas in regions of the Coulomb branch where there exists linear confinement we find 
that there exist stable string configurations for certain ranges of values of the parameter 
a . We finally discuss the case of open strings with fixed endpoints propagating in Rindler 
space, which turns out to have an exact classical-mechanical analog. 



Contents 



1 Introduction 2 

2 Wilson loops in AdS/CFT 4 

3 Stability analysis 6 

3.1 Small fluctuations 6 

3.2 Zero modes 8 

3.2.1 Transverse zero modes 9 

3.2.2 Longitudinal zero modes 10 

3.2.3 Angular zero modes 10 

4 Applications 13 

4.1 Boosted non-extremal D3-branes 13 

4.2 Spinning D3-branes 15 

4.2.1 Two equal nonzero angular momenta 16 

4.2.2 One nonzero angular momentum 18 

4.3 Marginally-deformed D3-branes 20 

4.3.1 The conformal case 20 

4.3.2 The sphere 23 

4.3.3 The disc 28 

5 Discussion and concluding remarks 30 
A Angular Schrodinger potentials 31 
B Angular fluctuations and the Lame equation 34 
C Soap films and strings in Rindler space 37 



1 



1 Introduction 



The AdS / CFT correspondence [1] maps the computation of the Wilson-loop heavy quark- 
antiquark potential in the planar limit of M = 4 SYM to the classical minimal-surface 
problem of calculating the action of a string connecting the quark and antiquark on 
the boundary of AdSs and extending into the radial direction [2]. However, when this 
method is applied to more general backgrounds [3, 4, 5, 6], one often encounters behaviors 
that are in sharp contrast with expectations based on the gauge-theory side [4]. A 
possible resolution could be that the parametric regions giving rise to these behaviors 
are perturbatively unstable and hence unphysical. 

Let us summarize the basic relevant facts for the case of M = 4 SYM. In the stan- 
dard conformal case, dual to a stack of D3-branes, this computation yields the expected 
Coulomb potential and all fluctuations about the classical string configuration are found 
to be stable [7, 8]. In extensions to the theory at finite temperature and at the Coulomb 
branch, dual to non-extremal and multicenter [9, 10, 11, 12] D3-branes respectively, one 
expects to find a screened Coulomb potential typical of thermal and Higgsed theories. 
However, one actually encounters situations where (i) the potential has a second branch 
of higher energy than the first, (ii) the potential exhibits a confining behavior at large 
distances, and (iii) the screening length is heavily dependent on the location of the probe 
string in the internal space. Although these types of behavior are quite counterintuitive, 
in a recent work [13] we established the reassuring result that the corresponding para- 
metric regions represent string configurations that are perturbatively unstable, with the 
physical regions giving indeed the expected screened Coulomb potential. Our analysis, 
partly motivated by a mechanical analog, was based on the zero-mode behavior of the 
differential equations governing the small fluctuations about equilibrium and led to a for- 
malism by means of which instabilities may be found by exact or approximate analytic 
methods. An important fact emerging from this analysis was the existence of instabilities 
due to fluctuations of non-cyclic angular coordinates, often appearing in cases where they 
are not a priori expected. In particular, these fluctuations cast the parametric region 
where confinement appears as unstable. 

The analysis of [13] was restricted to diagonal metrics, where all fluctuations obey de- 
coupled differential equations. However, there are physically interesting situations where 
the gravity duals where Wilson loops are evaluated contain non-diagonal metrics. A first 
class of such backgrounds are boosted (non-extremal) D3-brane backgrounds, used for 
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evaluating the potential for a quark- ant iquark pair moving with respect to the thermal 
medium [14, 15, 16, 17, 18, 19, 20]. For this case, the potential has a similar form to that 
in the zero-boost case, namely it is a double-branched function of the length with the 
lower branch corresponding to a screened Coulomb potential, but the screening length 
falls off as the velocity is increased. For this background, a numerical stability analysis 
has been presented in [21], indicating that the upper branch is perturbatively unstable. 
A second class of such backgrounds are spinning non-extremal D3-branes [9, 22, 23, 24], 
dual to J\f = 4 SYM at finite temperature and R-charge chemical potentials. The behav- 
ior is similar to that in the absence of chemical potentials, but the nontrivial dependence 
of the metrics on certain angles again calls for a stability analysis. Finally, a third class 
of such backgrounds are the Lunin-Maldacena deformations [25] (see also [26, 27]) of D3- 
brane backgrounds, dual to the Leigh-Strassler [28] deformations of M = 4 SYM which 
break supersymmetry down to M = 1. These backgrounds are characterized by the two 
real parameters 7 and a. In the usual situation where the quark and antiquark are not 
taken to be separated in the internal space, only the cx-part of the deformation affects the 
potential (see [29] for investigations on the effect of 7-deformations), resulting in various 
behaviors ranging from the standard Coulomb behavior (for deformed AdS 5 x S 5 ) to com- 
plete screening and linear confinement [30, 31] (for the deformed multicenter D3-branes). 
To investigate the significance of these results, a stability analysis is again in order. We 
note that now the appearance of a confining potential is not in conflict with expectations 
from the gauge-theory side. 

This article is organized as follows: In section 2, we give a brief review of the evaluation of 
Wilson loops in the backgrounds under consideration. In section 3, we analyze small fluc- 
tuations about the classical string configurations and we establish general results which 
allow us to determine the regions of instability using exact and approximate methods, 
generalizing the results of [13] to the non-diagonal case. In section 4, we apply these 
results to the gravity duals under consideration and we identify all unstable regions. In 
section 5, we summarize and conclude. In appendix A, we present the complete list of 
angular Schrodinger potentials for the backgrounds under consideration. In appendix B, 
we present the detailed solution of the equation for the angular fluctuations for a special 
case, in order to provide a consistency check of the numerical calculation employed in 
subsection 4.3.2. In appendix C, we apply our results to the problem of open strings in 
Rindler space, which turns out to have an exact classical-mechanical analog, namely the 
problem of the shape of a soap film stretched between two circular rings. 
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2 Wilson loops in AdS/CFT 



In the framework of the AdS/CFT correspondence, the calculation of the Wilson- loop 
potential of a heavy quark-antiquark pair proceeds by considering a fundamental string 
whose endpoints lie on the two temporal sides of the Wilson loop, and which extends in 
the radial direction of the dual supergravity background so as to extremize its worldsheet 
area. This type of calculation was first considered in [2] for the conformal case and was 
extended to more general cases in [3, 4, 5, 15, 29]. Below, we give a brief review of this 
procedure, adapted to the case where the metric has off-diagonal elements in the direc- 
tions transverse to the quark-antiquark axis and the 5-field has nonzero components. 

We consider a general background specified by a metric of the form 

ds 2 = G a dt 2 + 2G u dtdxi + Gijdxidxj + G yy dy 2 + G uu du 2 + G ab d9 a d9 b + . . . , (2.1) 

where, y denotes the (cyclic) coordinate along the spatial side of the Wilson loop, u 
denotes the radial direction extending from the UV at u — > oo down to the IR at some 
minimum value u min determined by the geometry, X; L stands for a generic cyclic coordinate, 
and 9 a stands for a generic non-cyclic coordinate. We also consider a 5-field of the form 

B 2 = B ai d6 a A dxi , (2.2) 

as is the case with many interesting gravity duals of gauge theories [25, 31, 32, 33, 34]. 
For future convenience, we introduce the functions 

g{U, da) = —G tt G UU , fy(U, 6 a) = —G tt Gyy , 

fij(u,9 a ) = GuGtj — G u Gij , fab(u,9 a ) = —G a G a b , h(u,9 a ) = G yy G uu . (2.3) 

According to AdS/CFT, the potential energy of the quark-antiquark pair is given by 

e~ iET = (W(C)) = e iS W , (2.4) 

where S[C] is the action for a string propagating in the supergravity background whose 
endpoints trace the contour C. The latter is given by the sum of Nambu-Goto and 
Wess-Zumino terms, 

S[C] = ~J drda (y-det^-ie^) , (2.5) 

where g a p and b a p stand for the pullbacks 

9 a p = G^d a x^dpx u , b a(3 = B^daX^dpx" . (2.6) 
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To proceed, we employ the gauge fixing 

t = r, u = a, (2.7) 
we assume translational invariance along t, and we consider the radial embedding 

y = y(u) , Xi = const. , 6 a = aO = const. , (2.8) 
supplemented by the boundary condition 

U (=4) = °° ' (2 ' 9) 
appropriate for a quark placed at y — —L/2 and an antiquark placed at y — L/2. In the 
ansatz (2.9), the constant values 9 a0 of the non-cyclic coordinates 9 a must be consistent 
with the corresponding equation of motion. As we shall see later on, this requires that 

d 

d a g(u, e a )\ da= e a0 = d a f y (u, 6 a )\ 9a =e a0 = ; d a = — . (2.10) 

For the above ansatz, only the Nambu-Goto part of the action is nonzero, 1 leading to 

J duy/g(u) + f y (u)y» , (2.11) 

where T denotes the temporal extent of the Wilson loop, the prime denotes a derivative 
with respect to u while g(u) = g(u,9 a0 ) and f y (u) = f y (u,6 a o)- Conservation of the 
momentum conjugate to y implies that the classical solution y c i{u) satisfies 

y' cl = ±^ , (2.12) 

where u is the value of u at the turning point, f y0 = f y (u ), the two signs correspond 
to the two branches around the turning point, and F is defined as 

F = 9fy f . (2.13) 

Jy - JyO 

Integrating (2.13), we express the separation length as 

POO a/P 

L = 2f l y ' 2 du^-, (2.14) 
Juo Jy 

while inserting (2.12) into (2.11) and using (2.4), we obtain the potential energy 

E = - duVF-- du^fg . (2.15) 

When the integrals (2.14) and (2.15) can be evaluated exactly and the first one can be 

inverted for uq, these equations lead to an explicit expression for E = E(L). However, in 

practice this cannot be done, except for a few simple cases, and Eqs. (2.14) and (2.15) 

are rather regarded as parametric equations for L and E with parameter uq. 

1 Nonzero contributions from the Wess-Zumino part can arise when B tu and/or B ty are nonzero (see 
[6]). However, this occurs in a very restricted class of backgrounds. 
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3 Stability analysis 



We now turn to a stability analysis of these configurations just discussed, our goal being 
to identify all parametric regions which are unstable and for which information obtained 
using the gauge/gravity correspondence cannot be trusted. In this section we generalize 
the results of [13] to non-diagonal metrics of the form (2.1) and we establish a series of 
results which will ultimately allow us to identify the stable and unstable regions by a 
combination of exact and approximate analytic methods. 



3.1 Small fluctuations 

To investigate the stability of the string configurations of interest, we consider small 
fluctuations about the classical solutions discussed above. In particular, we will be in- 
terested in three types of fluctuations, namely (i) "transverse" fluctuations, referring to 
the cyclic coordinates Xi transverse to the quark- ant iquark axis, (ii) "longitudinal" fluc- 
tuations, referring to the cyclic coordinate y along the quark- ant iquark axis, and (iii) 
"angular" fluctuations, referring to the non-cyclic coordinates 9 a . The above fluctuations 
may be parametrized by keeping the gauge choice (2.7) unperturbed and perturbing the 
embedding as 

Xi = 6xi(t,u), y = y d (u) + 6y(t,u) , 9 a = 9 a0 + 69 a (t, u) . (3.1) 

Inserting this ansatz in the action (2.5) and expanding in powers of the fluctuations, we 
obtain the series 

S = S + S 1 + S 2 + . . . , (3.2) 

with the subscripts corresponding to the respective powers of the fluctuations. The 
zeroth-order term gives just the classical action. The first-order contribution reads 



Si = 

2tt 



J dtdu y/fo by + J ^GuSxi + {^^d a g + ^°^ 2 d a f^j 59 a 



, (3.3) 



where d a g = d a g(u,9 a )\e a= Q a0 and d a f y = d a f y (u,9 a )\o a= o a0 . The first term is a surface 
contribution which is exactly cancelled by a similar term with the opposite sign corre- 
sponding to the contribution of the linear order action coming from the fluctuations of 
the lower string (cf. (2.12)), provided we keep the variations of both strings at u — Uq 
equal. The second term is a total time derivative and hence is completely irrelevant. 
Finally, in the third term, the coefficient of 59 a is just the equation of motion of 9 a and 
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the requirement that it vanish leads indeed to the conditions (2.10). The second-order 
contribution is written as 



2tt 



dtdu 



2FV2 



2j p3/2 



2j pl/2 



9fy 

( 1 



G ti (5y'5xi - fo)<^) + B ai (59'Jxi - <W„&r-) 



(3.4) 



a 6 



hF 1 ' 2 
2gf y 



se a 59 b 



+ 



' 9 a d b g + f -y^d a d b f y ) 60.M 



4 4 2 



again with all functions and their 6* a -derivatives evaluated at 9 a = 9 aQ . Writing down the 
equations of motion of the fluctuations and introducing a harmonic time dependence, 



5x»(t,u) = 5x»(u)e 



-iu)t 



(3.5) 



we obtain the equations 



du \F l / 2 du 
d ( gf y d 



d ( fij d \ 1 



gfy 



5Xj — ILO 



d u \G u J^)5y-d u B ai 59 a 



, 



du \F 3 / 2 du 

d f fab d 



h 



du \ F 1 / 2 du 



F 1 / 2 
2 hF 1 / 2 f l 



5y + iud u G t 



gfy 



5xi = , 



(3.6) 



d a d b g + %^ d a dbfy 



2FV2 



59 b - \ujd u B ai 5xi = 



gfy v-- -jy 

for the transverse, longitudinal and angular fluctuations respectively. We see that the 
fluctuations generically couple to each other, satisfying a system of equations of the form 



d_ 

du 



f p x \ 

Py 

V o op,/ 



d 

du 



( u 2 q x -iuj~f 







Sy 

\se J 



0, (3.7) 



\ -iu}(3 u 2 q e + r e ) 

where the matrices p^ = (p^-), p e = (p ab ), q x = (<^), q e = (q ab ), r e = (r ab ), (3 = (f3 ai ) 
and the column vector 7 = (7^) are read off from Eqs. (3.6). The problem is defined in 
the interval 

u ^ u < 00 , m > u min , (3.8) 
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and the boundary conditions for the fluctuations, determined by requiring that the string 
endpoints be held fixed and that the two parts of the string glue smoothly at u = Uq read 

lim 5xi(u) = , lim (u — uqY^Sx'^u) = , 

lim Sy(u) = , lim [8y(u) + 2(u - u Q )5y\u)\ = , (3.9) 

lim 56 a {u) = , lim (u - u ) 1/2 56' a (u) = . 

These boundary conditions follow from a straightforward extension of the proof given in 
[13] and is essentially based on the fact that u = w and u = oo present regular singular 
points of the system (3.6). Therefore, our stability analysis has reduced to a coupled 
system of Sturm-Liouville equations, with our objective being to determine the range of 
values of u for which the lowest eigenvalue becomes negative. 

Ideally, one would like to solve the above Sturm-Liouville equations exactly, obtain the 
lowest eigenvalue Uq in terms of the parameter u and determine the regions where Uq 
becomes negative. However, in most cases, this is impossible due to the complexity of 
the equations. On the other hand, it turns out that we can obtain useful information by 
studying a simpler problem, namely the zero-mode problem of the associated differential 
operators. Regarding the transverse and longitudinal fluctuations, we will use our Sturm- 
Liouville description to prove that transverse zero modes do not exist while longitudinal 
zero modes are in one-to-one correspondence with the critical points of the function L(u ). 
For the angular fluctuations, we will employ an alternative Schrodinger description which 
allows us to identify angular zero modes using either exact or approximate methods. 



3.2 Zero modes 

To motivate the significance of zero modes for our stability analysis, we consider a Sturm- 
Liouville system of the form considered earlier on and we let uj^(u ) be the corresponding 
eigenvalues. By standard results of Sturm-Liouville theory, these eigenvalues are real 
and strictly-ordered in the sense that cc;q(-u ) < ^i( M o) < • • •, from which it also follows 
that different eigenvalues do not cross as we vary u . Furthermore, we know that, for 
sufficiently large values of u , all eigenvalues are positive. The above considerations imply 
that the first occurrence of instabilities will arise at a value of u at which the lowest 
eigenvalue Uq(uo) becomes negative. Therefore, a necessary condition for the appearance 
of instabilities is the existence of values of u for which Wq^o) — 0. To verify that 
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these zero modes really correspond to points where uJq changes sign, we may solve the full 
Sturm-Liouville system near each critical point v,®, a task that may be accomplished 
using perturbative [13] or numerical methods. In fact, as we shall see later on, in all 
examples considered in this paper, a zero mode always signifies a change of sign of uIq 
i.e. marks the boundary between a stable and an unstable region. 

Restricting to zero modes simplifies our problem tremendously, as the three equations in 
(3.6) decouple, reducing to 

d_ ( fij d 
du \F l l 2 du" 3 



d ( gf y d 



du \ F 3 / 2 du 



5xj J = , 

Sy)=0, (3.10) 

56 b = . 



Tu " {^ dadb9 + ^p- dadbfy ) 



V 2 y 

Using these simplified expressions, we will next prove that transverse zero modes do not 
exist and that longitudinal zero modes are in one-to-one correspondence with the critical 
points of the function L(u ), and we will devise analytic methods for seeking angular 
zero modes. This is a fairly straightforward extension of the results of [13]. 

3.2.1 Transverse zero modes 

We consider first the case of the transverse fluctuations. Using the definition of F in 
(2.13), performing an integration by parts, and expanding about u = uq, we write the 
general solution of the first of (3.10) as 

5xi = [ duj 9 ^ v f^Cj 

Ju y Jy ~ JyO 

= ■ v'TWW:, ; '•, _/ ; ' X <>" v'TWVA, [^jhl.j 1 j '•, (3.H) 

= -2 r du^f-^o d u (^f-A c, - 2. f^f-^iu - Mo ) 1/2 + 0(u - u ) , 

Juo \ Jy J y JyO 

where the q are constants. In order for this zero mode to exist, it must satisfy the first 
boundary condition in (3.9), which requires the coefficient of (u — Uq) 1 ^ 2 to vanish. As the 
matrix f^l turns out not to admit any nontrivial null eigenvectors, 2 this is only achieved 



2 We have not a general proof of that statement, but this is indeed the case in all examples of the 
present paper as well as of [13]. 
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when all q = i.e. for the trivial solution 5xi = 0. Therefore, the transverse fluctuations 
do not lead to zero modes of the system (3.6). 



3.2.2 Longitudinal zero modes 

Turning to the longitudinal fluctuations, the fact that they decouple from the rest in the 
zero-mode analysis implies that the results of [13] for the diagonal case hold in our case 
as well. In particular, in [13] it was proven that the longitudinal zero-mode solution can 
be expressed in terms of the wo~derivative of the length function as follows 



*v ~ 2 \Hf ( M - u ^ 1/2 + + o ((« - Mo ) 1/2 ) • (3.i2) 

V Jyo ho 

In order for this zero mode to exist, it must satisfy the second boundary condition in 
(3.9), which requires the constant term to vanish. Therefore, the solution exists only if 

L'(u Q ) = , (3.13) 

or, equivalently, if [13] 

*• d, 1=0, (3.14) 



1 \f7~u .///ii y fy 

which requires that the derivative term should change sign at least once as u ranges in 
the interval u e [uq,oo). Therefore, given the critical points of the length function, we 
may determine all values uq c of uq where the system (3.6) has a zero mode. In view 
of our earlier remark that zero modes always signify a transition from a stable to an 
unstable region in our examples, the longitudinal instabilities for the various configura- 
tions of interest are as shown in Fig. 1. Additional instabilities may arise from angular 
perturbations, so that part(s) of the curves, stable under longitudinal perturbations only, 
might be actually unstable. 



3.2.3 Angular zero modes 

We finally consider angular fluctuations, which must be treated separately due to the 
fact that the non-cyclicity of the angular coordinates induces a "restoring force" term 
in the corresponding Sturm-Liouville equation, implying that we cannot write 56 a in an 
explicit integral form, even in the zero-mode case. On the other hand, if we restrict to 
situations where the angular fluctuations are decoupled from the rest (which is indeed 
the case in most of our examples), the equation satisfied by each of the fluctuations can 
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(a) 



(b) 



(c) 



(d) 



Figure 1: The four generic types of behavior of L(u ) and E(L) encountered in AdS/CFT 
calculations of Wilson loops, corresponding to (a) no critical points, Coulomb or screened 
Coulomb potential, (b) no critical points, Coulomb/confining potential, (c) one critical 
point, double-branched potential and (d) two critical points, multi-branched potential. 
The dashed parts of the curves indicate regions that are unstable under longitudinal 
perturbations. 

be easily brought to a Schrodinger form. This description allows us to determine the zero 
modes to quite high accuracy by using approximate methods. Following [13], we outline 
the procedure below. 

Letting 6 be any of the decoupled non-cyclic angular variables, the Sturm-Liouville equa- 
tion for its fluctuations reads 




(3.15) 



Employing the change of variables 





(3.16) 



we may transform Eq. (3.15) to a standard Schrodinger equation 




(3.17) 



with the potential 




+ (poqe)- 1/4 -f,( P oqe) 1/4 ■ (3.18) 
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The problem is defined in the interval 

roc r^- 

^ x ^ x , x = / duA — , (3.19) 

J no \ Pe 

which is finite, making the fact that the fluctuation spectrum is discrete manifest. The 
boundary conditions are simply 

tf(0) = , *'(xo) = . (3.20) 

Note that, in general, Eq. (3.16) does not lead to a closed expression for u in terms of x 
and therefore it is not always possible to write down the potential as an explicit function 
of x. In what follows, we will encounter cases where Vg(x) can be determined exactly as 
well as cases where it may be adequately approximated by solvable potentials. 

To develop our approximate methods, we first note that the behavior of the angular 
Schrodinger potential Vg(u; u ) in the limits u — > oo and u = u is given by 

V e (oo,u ) =au 2 + V OQ , (3.21) 

and 

V e (u ; u ) = - / 9 U0 — ^ + -<% = ^0 , (3.22) 

where a, and V are finite quantities that depend on the parameter u . When a ^ 0, 
the potential expressed in terms of the variable x rises from a minimum value to infinity 
in the finite interval x G [0, x ], and hence it is reasonable to approximate it by an infinite 
well given by 

K„™ = ( 7 - ° *' ) ; vJ ^ + a = \. (3.23) 

I 00 , otherwise j j V , a^O j 

With the boundary conditions (3.18) the energy levels read 

4M« )] 2 

When a = 0, this approximation is valid for both the ground state and the excited states. 
When a > however, the approximation is valid only for the ground state as excited 
states are affected by the details of the exact potential. In any case, we may determine the 
critical value -u c by solving the equation ^(^o) = 0; clearly, a solution to this equation 
exists only if the average V is negative at least in a finite range of values of uq. The 
infinite-well approximation just described is to be used for obtaining an indication for 
the existence of zero modes and for estimating u 0c . In most cases, these estimates are 
quite close to the exact values, determined through numerical analysis. 
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4 Applications 



Having developed the necessary formalism, we now turn to analyzing the stability prop- 
erties of the string configurations of interest in backgrounds of boosted non-extremal 
D3-branes, spinning D3-branes and marginally-deformed D3-branes, all of which fall in 
the general category described by the metrics (2.1). For each case under consideration, 
we give the formulas for the length and the potential energy and we perform a stability 
analysis according to the guidelines of the previous section. Depending on the problem 
at hand, the regions of stability are determined by exact or approximate analytic meth- 
ods, by the numerical solution of certain algebraic or transcendental equations, or by 
the numerical evaluation of certain integrals. In all cases, this represents a considerable 
improvement, both conceptual as well as practical, over the direct numerical solution of 
the differential equations governing the fluctuations. 

4.1 Boosted non-extremal D3-branes 

The first type of metrics we will consider are obtained by a boost of the metric for non- 
extremal D3-branes along one brane direction transverse to the quark-antiquark axis, say 
x. They are given by 

_ (l _ f£) ^ + 2 1 ^dtd X +(i + t^] dx * + dy i + d2 ? 
+ R 2 l-^-i du 2 + dnl) , (4.i) 

where v is the boost velocity and 7 = 1/ yl — v 2 . The metric has a horizon at 

u H = n, (4.2) 

while there also exists a velocity-dependent radius w 7 = y/jfJ, ^ u H beyond which the 
string cannot penetrate (see [14, 15, 16, 17, 19] for discussions). The Hawking temper- 
ature of the solution is given by T = [i/itR 2 . For the stability analysis for this metric, 
the "transverse" coordinates are (x, z, f2 5 ), with x coupled to the longitudinal coordinate 
y, while no "angular" coordinates exist. This latter fact implies that the position of the 
string in the internal space can be arbitrary. 

The calculation of the Wilson loops of section 2 for the above metric yields the potential 
for a quark and antiquark moving with velocity v with respect to a thermal plasma at a 



u 



13 



1.5 
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Figure 2: Stability diagram in the (v, Uq) plane for the case of the boosted non-extremal 
D3-brane background. The light-shaded region corresponds to instabilities under longi- 
tudinal perturbations and the dotted curve indicates our approximation (4.7) for woc(7)- 
The dark-shaded region corresponds to the values u < ^7 which are inaccessible to the 
string. The region u < 1 beyond the horizon has been excluded from the diagram. 

temperature T and may be used as a crude model for the study of meson dissociation in 
plasmas. Switching to dimensionless units by the change of variables 



u 



L — > —L 



we find that the length and the potential energy read [15, 19] 
and 



2y/27r 3 / 2 y/?4 - 7 2 

r(i/4) 2 ui 



1 1 2'4'4'ug 



y^TT 3 / 2 

r(i/4) 2 



«0 2^1 



1111 

"I' 2' Z'ut 
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4' 


4' 





+ 1 



(4.3) 



(4.4) 



(4.5) 



where 2-^1(0, b, c; x) is the standard hypergeometric function. The behavior of the length 
and the energy is as in Fig. lc, with the maximal value of the length, L c (1x0,7), being 
a decreasing function of the velocity and satisfying the approximate law L c (u ,7) — 
7~ 1 / 2 L c (-u , 1) [15] indicating an enhancement of the dissociation rate of the quark- 
antiquark bound state with increasing velocity. 

Turning to the stability analysis, our general results imply that the coupled system of 
the (Sy,Sx) fluctuations has instabilities for uq < Moc(7) where woc(7) is the critical 
point of L(uq, 7)- The latter is determined by the solution of (3.13) which leads to the 
transcendental equation 



5u 4 (u 4 - 37 2 ) 2^1 



13 5 1 



6(4-^)2^ 



3 7 9 1 

2'4'4'^f 



0, 



(4.6) 
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which can be obtained numerically. For the physically interesting case 7 ^> 1 (in which 
we also have u ^> 1) we can solve this equation perturbatively in I/7 with the result 

1 23 



:(7) = 3 1/4 v/7 



0(7" 



(4.7) 



157 2 13507 4 

The results of the stability analysis just presented are summarized in the diagram of 
Fig. 2. Note that the estimate (4.7) for wo c (7) is remarkably close to the exact numerical 
result even for low velocities. Our results are in accordance with those of the numerical 
analysis presented in [21]. 

The above analysis may be extended to the case where the velocity and the axis of 
the quark- ant iquark pair are not perpendicular but form an angle smaller than 7r/2. 
For such a case, we may modify our ansatz (2.8) to include a variable x that has non- 
trivial dependence on u, while in the special case of motion parallel to the axis we 
may equivalently boost along the y instead of the x axis which would mean that the 
metric (2.1) would include an extra G ty term. Such string solutions were considered in 
[15, 17, 19, 20] and the behavior of the length and energy is again as in Fig. lc. Hence 
we expect instabilities under longitudinal perturbations for -u < -u 0c (7), a fact that we 
have actually verified by a small-fluctuation analysis which is however too lengthy to be 
included here. 



4.2 Spinning D3-branes 

We next consider the case of spinning (non-extremal) D3-branes, dual to M = 4 SYM 
theory at finite temperature and R-charge chemical potentials. These metrics were found 
in full generality in [9], based on previous results from [22], and their thermodynamical 
properties were examined in [24]. In the conventions of [23] which we here follow, the 
field-theory limit of these solutions is characterized by the non-extremality parameter 
/i and the angular momentum parameters Oj, i = 1, 2, 3. Here, we will restrict to two 
special cases, corresponding to two equal nonzero angular momenta, a 2 = a 3 = r , and 
one nonzero angular momentum, ai = r , to which we will apply our stability analysis. 
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4.2.1 Two equal nonzero angular momenta 



For the case of two equal angular momenta, the metric reads 



ds 2 = H- 1 ' 2 
+ H 1 ' 2 



- 1- 



dt z + dx z + dy z + dz z 



[U 4 — /i 4 )(w 2 — Tq) 



[u 2 - r 2 cos 2 0)d6 2 + u 2 cos 2 0dQ 2 3 + {u 2 - r 2 ) sin 2 #d0 2 



(4.8) 



' i? 2 



(it cos 2 6*(sin 2 -?/><i02 + cos 2 ipdfe) 



where 



H 



R 4 



u 2 {u 2 — r 2 cos 2 0) 

and dill = dif) 2 + sin 2 ipd§\ + cos 2 -0^03 • The metric has a horizon at 

u H = H, 



(4.9) 



(4.10) 



and considerations of thermodynamic stability restrict the ratio A = /i/r according to 

A ^ 1 (CE) , A ^ V2 (GCE) , (4.11) 

where the two values refer to the canonical and grand-canonical ensembles respectively, 
with the former giving the lowest bounds in our examples. 3 For the stability analysis 
for this metric, the "transverse" coordinates are (x, z, ijj, $1,2,3), with 2 ,3 coupled to 
the longitudinal coordinate y, and the only "angular" coordinate is 0. The restriction 
(2.10) on the angular location O of the string allows only trajectories with O = and 
#0 = 7r/2. To keep the discussion to a reasonable length, we will here restrict to the 
6> = vr/2 trajectories. For the trajectories with O = 0, the longitudinal and angular 
fluctuations may be analyzed in the same way, using the Schrodinger potentials given in 
appendix A for the latter case. 

To examine the Wilson-loop computation, we switch to dimensionless units by setting 



u — > r u , 



u -> r u , 



E — > —E . 

IT 



(4.12) 



Then, we find that the length and the potential energy are given by the expressions 

duu 



V /( M 2_ 1)(M 4_ M 4 )(M 4_ A 4) ' 



(4.13) 



3 For details on the thermodynamic stability of spinning branes see [24]. For a brief summary for the 
type of the spinning D3-branes used in this paper see [35] . 
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1.14 




Figure 3: Stability diagram in the (A, Uq) plane for the spinning D3-brane with two angu- 
lar momenta and for the trajectory with 6q = n/2. The two light-shaded regions on the 
upper right-hand and lower left-hand sides correspond to instabilities under longitudinal 
and angular perturbations, respectively. The dark-shaded region corresponds to values 
of u beyond the horizon. The region to the left of the vertical dashed line corresponds to 
the values A < 1 for which the background is thermodynamically unstable. The dotted 
line corresponds to the curve drawn using the approximation described below (4.15). 



and 



A 4 



u 



l 



(4.14) 



which unfortunately cannot be evaluated in closed form. 4 Their behavior is as in Fig. lc. 

Turning to the stability analysis, our general results imply that the coupled system of 
the (<5?/, 502,3) fluctuations has instabilities for u < u^J(X) where Mq^A) is the critical 
point of L(u , A). This value can be obtained by numerically solving the equation (3.14) 
which after some algebra takes the simple form 

u 6 - 3X 4 u 2 + 2A 4 



du- 



. 



(4.15) 



UQ 



«V(« 2 -i) 3 (« 4 -^)(« 4 -a 4 ) 

In order for a solution to exist, the integrand must change sign at least once in the 



interval u G [uq,oo). This requirement leads to the inequality A < uq < A a/2 cos(0/3), 
where cos0 = — 1/A 2 , hence G [7i/2,n]. Taking the upper bound as an equality, i.e. 
u^J(X) = A a/2 cos(0/3), gives a quite good approximation to the exact numerical result, 
as one can see by inspection of Fig. 3. 



4 Such integrals can be thought of as periods of Riemann surfaces corresponding to algebraic curves. 
In this case the genus of the Riemann surfaces is at least two. 
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Finally, for the angular fluctuations, the Schrodinger potential reads 



V e = -l, 



(4.16) 



implying that the infinite-well approximation is exact in this case. Therefore, the critical 
value of -uo beyond which instabilities occur are obtained by using Eq. 3.24, which is in 
fact exact in this case, with V — — 1 and solving the equation to 2 = 0. The resulting 
equation has the form 

duu 3 



f 

Jun 



7T 



(4.17) 



, U0 vV-l)K-«o)K-A 4 ) 2 ' 

and again can be solved numerically to give w^A). The results of this stability analysis 
are summarized in the diagram of Fig. 3. 



4.2.2 One nonzero angular momentum 



For the case of one angular momentum, the metric reads 



ds 2 = H- 1 ' 2 
+ H 1 ' 2 



fi 4 H 



dt 2 + dx 2 + dy 2 + dz 1 



+ 



ffl/2 u 2 (u 2 + r 2 cos 2 6) ^ 2 



u 4 + r 2 u 2 — /i 4 



(u 2 + r 2 cos 2 6)d6 2 + u cos 2 6dn 2 3 + (u 2 + r 2 Q ) sin 2 6d(f)i 



(4.18) 



where 



- 2^ sin 2 edtdfa 



H 



R 4 



u 2 (u 2 + Tq cos 2 6) 



and dVt\ is as before. This metric has a horizon at 

u 2 n = \ {-r 2 + V^o 4 + 4/i 4 ) , 
and thermodynamic stability restricts A = /i/r to the range 

A > 0.685 (CE) , A > 0.93 (GCE) , 



(4.19) 



(4.20) 



(4.21) 



for the canonical and grand-canonical ensembles respectively. For the stability analysis, 
the "transverse" and "angular" coordinates are as before, but now it is 0i that couples 
to the longitudinal coordinate y. Again, the allowed trajectories have 9 = or 9 Q = ir/2, 
and only the latter case will be considered here. 
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0.69 1 



Figure 4: Stability diagram in the (A, Uq) plane for the spinning D3-brane with one angu- 
lar momentum and for the trajectory with 0q = it/ 2. The light-shaded region corresponds 
to instabilities under longitudinal perturbations. The dark-shaded region corresponds to 
values of uq beyond the horizon. The region to the left of the vertical dashed line corre- 
sponds to the values A < 0.685 for which the background is thermodynamically unstable. 
The dotted line corresponds to the curve drawn using the approximation described below 
(4.24). 



Using the same rescalings as before, we find that the expressions for the length and energy 
read 

L(u , A) 



2\lul - A 4 



du 



\J (w 4 — Mq)(m 4 + m 2 — A 4 ) ' 



(4.22) 



and 



E(u ,X) 



du- 



u 4 -X 4 



du 



u ' 



A 4 



(4.23) 



„„ vV-^K + ^-A 4 ) Ju H V« 4 + « 2 -A 4 
which again cannot be evaluated in closed form. Their behavior is as in Fig. lc. 

Turning to the stability analysis, the coupled system of the (Sy, 5<pi) fluctuations has 
instabilities for for u Q < -Uoc(A) where Moc(A) is the critical point of L(uo, A), obtained by 
numerically solving Eq. (3.14) which for our case reads 

u 8 - 4A 4 m 4 - 4AV + 3A 8 



du- 







(4.24) 



, uo uV(u 4 -u 4 )(u 4 + u 2 -A 4 ) 3 
This equation has a solution only if the integrand changes sign at least once in the 
interval u G [u ,oo). This leads to the inequality uu < u < u nmx (X) where ■u max (A) is 
the largest root of the numerator for which we will not present an explicit expression 
due to its complexity. As before, saturation of the upper bound gives a curve defined by 
M oc(A) = w max (A) which gives a quite good approximation to the exact numerical result, 
as one can see by inspection of Fig. 4. 
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Finally, for the angular fluctuations, the Schrodinger potential reads 

V = l, (4.25) 

implying that no angular instabilities occur. The results of the stability analysis are 
summarized in the diagram of Fig. 4, which is qualitatively similar to Fig. 3 apart from 
the absence of angular instabilities. 

4.3 Marginally-deformed D3-branes 

We finally consider the Lunin-Maldacena deformations [25] of D3-brane solutions, dual 
to the Leigh-Strassler [28] marginal deformations of Af = 4 SYM which break supersym- 
metry down to J\f = 1. These backgrounds are characterized by the complex parameter 
(3 = j+ra (where r is the Type IIB axion-dilaton and 7 and a are real parameters) which, 
on the gauge-theory side, represents the complex phase entering the Leigh-Strassler su- 
perpotential (where r is the complexified gauge coupling). In particular, we will consider 
marginal deformations of the conformal AdSs x S 5 background [25] and of the multicenter 
backgrounds corresponding to D3-branes distributed on a sphere and on a disc [31]. 

4.3.1 The conformal case 

Starting from the deformation of the conformal AdSs x S 5 background, we rescale the 
deformation parameters as {j3, 7, a) — > jp(/3, 7, g s (T ) an d we write the deformed metric as 

ds 2 = H 1/2 



(4.26) 



^ + dx 2 + dy 2 + ^ +R 2(^_ + dQ 2 ^ 

where dill p is the metric on the deformed S 5 , given by 

dttlp = d6 2 + g sin 2 Odcft + cos 2 B\d^ 2 + £(sin 2 ^d<g + cos 2 ^#1)] 

+ g\(3\ 2 cos 4 6 sin 2 6 sin 2 2iP {dfa + d<f) 2 + d<f) 3 ) 2 , (4.27) 

and the functions g and 7i are given by 

g- 1 = 1 + 4\[3\ 2 cos 2 ^(sin 2 6 + cos 2 9 cos 2 tp sin 2 V) , 

U = 1 +4(j 2 cos 2 #(sin 2 # + cos 2 #cos 2 ^sin 2 ^) . (4.28) 

The solution also includes a nonzero -B-field which, in the case when the deformation 
parameter 7 vanishes, is equal to 

B 2 = - cos 4 9 sin 2ip(d<f>i + d<p 2 + d(f) 3 ) A dip , (4.29) 
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and hence is of the form (2.2). For the stability analysis for this and all subsequent 
metrics, the "transverse" coordinates are (x, z, ^1,2,3), with all of them being decou- 
pled from the longitudinal coordinate y but with the <pi coupled to each other, while 
the "angular" coordinates are (9,ip), and are decoupled from each other. The restric- 
tion (2.10) on the angular location (9o,ipo) of the string allows only for the trajectories 
(#0,^0) = (0,tt/4), (^cV'o) = (0,0 or n/2) (with both choices for tp leading to equiv- 
alent results), (9 ,^ ) = (tt/2, any), and (9 ,^ ) = (sin~ 1 (l/y / 3), vr/4). All trajectories 
but the last will also be valid for the multicenter case with the branes distributed on the 
sphere and on a disc. We also note that for the trajectories with 9 = n/2 the variable ip 
becomes a completely decoupled "transverse" coordinate and hence its fluctuations are 
stable. Also, the presence of the nonzero 5-field in (4.29) implies that, for the trajec- 
tories with (9o,ipo) = (0, 7r/4) and ^0) = (sin _1 (l/v^3), tt/4), the angular fluctuation 
Sip and the transverse fluctuations Sfa are actually coupled. However, the coupled terms 
are proportional to the eigenvalue u> and hence do not affect at all the results of our 
analysis which is based on the zero modes. For (9 0: ip ) = (0,0 or 7r/2), (7r/2,any), Sip 
and 5(f>i are decoupled. 

For the present case, the potential energy E can be calculated explicitly as a function of 
L with the result 



where k(a) is an angle- and cr-dependent factor, given explicitly in the examples below, 
that becomes unity for a = 0. The factor in parentheses is the result of [2], giving the 
standard Coulomb behavior expected by conformal invariance. We note in passing that 
this factor is unaffected by turning on 7-deformations, unless the quarks are given a 
separation in the internal deformed sphere as well [29]. 

Regarding stability, the transverse/longitudinal fluctuations are obviously stable, while 
the Schrodinger potentials for the angular fluctuations all have the form 



where ag^(a) are angle- and cr-dependent factors that can be read off the formulas of 
appendix A. The expression for the Schrodinger variable x in terms of u is 




(4.30) 



Ve^ = a e ^(a)u . 



(4.31) 




(4.32) 



and, accordingly, the value of the endpoint x is 

r(i/4) 2 1 

X = - nz — 

4V27T u 



(4.33) 
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Instabilities may occur only when one of ag^(a) is negative. Although Eq. (4.32) cannot 
be inverted in terms of u to allow for an analytic solution of the Schrodinger problem, 
we may obtain a lower bound for o above which instabilities definitely occur. We first 
note that the infinite-well approximation is not valid here, as the potential satisfies Eq. 
(3.21) with a < 0. To obtain our bound, we consider the limit u/uo 1, where x ~ 1/u 
and Ve^ — ag^(a)/x 2 . Then, standard arguments from quantum mechanics [36] show 
that for, ae^{a) < — |, this potential supports an infinite tower of negative-energy states 
and the solution becomes unstable for all u . This gives the desired bound on a. The 
results for the allowed trajectories are as follows: 

• (9 ,ip ) = (0, 7r/4). For this case we have 

2a 2 . . Aa 2 

iT^ ' Ma) = ~TTa 2 



k(a) = Vl + a 2 , a e {a) = — — - , a^(a) = -— — . (4.34) 



Instabilities occur only from the Sip fluctuations for a > 1/V15 ~ 0.258. Since in the 
UV all marginally deformed backgrounds approach (4.26), (4.27), the above bound is 
universal as long as the corresponding trajectory remains valid. This is indeed the case 
for the sphere and disc brane distributions we consider. 

• (9o,ip ) = (0,0 or 7r/2). Here, we have 

k(a) = 1 , ae(a) = a^(a) = 4a 2 , (4.35) 

and no angular instabilities occur. 

• (9 ,ipo) = (7r/2,any). Here, we have 

k(a) = 1 , a e (a) = Aa 2 , (4.36) 

and again no angular instabilities occur. 

• (6* ,-?/>o) = (sin~ 1 (l/y / 3), 7r/4). Here, we have 



Aa 2 8a 2 
k(a) = \/ 1 + — , a e {a) = a^a) = ~j^~ 2 ■ (4.37) 



Instabilities occur from both 59 and Sip fluctuations for a > a/3/28 ~ 0.327. This bound 
does not survive in the multicenter cases we consider below since the corresponding 
trajectory is no longer valid. 

Note that the existence of un upper bound for the deformation parameter a beyond 
which stability breaks down is reminiscent of an analogous fact for giant gravitons on the 
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pp-wave limit of the deformed background (4.26), (4.27), which exist only for values of 
a < 1/ y/l2. These giant graviton solutions can be thought of as integrated perturbations 
around their infinitesimal size. This implies that every order in the small size perturbative 
expansion of the giant graviton solution is unstable if the above bound is not respected. 
The different upper value for a in this case should be attributed to the fact that the 
probes in the giant graviton case refer to D3-branes and not to strings. 

4.3.2 The sphere 

We next consider the deformation of the background corresponding to D3-branes dis- 
tributed on a sphere of radius r . Switching to dimensionless units by using Eq. (4.12), 
we write the deformed metric as [31] 

ds 2 = H 1/2 \h^ 2 {-dt 2 + dx 2 + dy 2 + dz 2 ) + H^ 2 u2 '^°f 6 du 
y u r 

+ H 1 ' 2 [(u 2 - r 2 cos 2 6)d6 2 + Q(u 2 - r 2 ) sin 2 6d<f> 2 ] 

+ H 1/2 u 2 cos 2 6 [di) 2 + £(sin 2 i)d(& + cos 2 ^#§)] 

1/2 |/j| 2 ^ 2 ( M 2 -r 2 )cos 4 gsin 2 gsin 2 2^ 2 

+ H M 2_ r 2 cos 2fl (#1 + + d< W } '( 4 - 38 ) 



with 



r -l , , A \a\2 2g ( U ~ r l) sin g + U COS 6 COS t Sil1 t 

y — J- + 4 \p\ cos a 9 2 9 /1 ' 

u l — Tq cos 2 

2 2 (m 2 - r 2 ) sin 2 g + m 2 cos 2 cos 2 ^ sin 2 y> 
7t = l + 4(j cos - 5 > (4.39) 



and H is given by (4.9). We also note that here there is a nonzero 5-field which for 

7 = is proportional to the expression given in (4.29). Now, the restriction (2.10) on the 

(0 O , ^o) of the string allows only for the first three trajectories considered earlier, namely 

(0 O , ^o) = (0,vr/4), (0 O , ^o) = (0,0 or vr/2) and (0 o ,Vo) = (vr/2,any). 

We next examine the potentials arising in each case in turn: 

• (0o,V'o) = (0,vr/4). The length and potential energy read (see sec. 6.2 of [31]) 

L(u , a) = - = 2U ° [n(a 2 , k) - K(k)] , (4.40) 

- i+^)(2u ~ T+js) 



Uf 



and 

E(u , a) = VTT^ | ^2ul--^— 2 [a 2 K(k) - E(fc)] + E(c) - Y^K(c) |> , (4.41) 
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where 



k = 




l+a 2 



1+CT 2 



1 



a = 




1 



(4.42) 



Note that these results are independent of the deformation parameter 7. The reason 
is that, as remarked also in the conformal case, we have not separated the quarks in 
the internal deformed-sphere space. Wilson-loop potentials where such a separation was 
considered, but with vanishing cr-deformation, were computed in [29]. 

For a — 0, the behavior of the length and energy is the same as for the undeformed 
case (Fig. lc). As a is turned on, the length and energy curves are reminiscent of van 
der Waals isotherms for a statistical system with uq, L and E corresponding to volume, 
pressure and Gibbs potential respectively (see e.g. [37]). In particular, there exists a 
critical value of a, given by ~ 0.209 (found analytically in [31]), below which the 
system behaves like the statistical system at T < T cr (Fig. Id) and above which the 
system behaves like the statistical system at T > T cr (Fig. lb). For nonzero a, we have 
the usual Coulombic behavior for large uq (small L), while in the opposite limit, uq — > 1, 
(large L), the asymptotics of (4.40) and (4.41) lead to the linear potential 



To examine the significance of these results, it is crucial to examine the stability of the 
corresponding string configurations. 

Our general results imply that the solution is stable under transverse perturbations. For 
the longitudinal fluctuations, the fact that L(u ) has two extrema u^\a) and u^\a) 
for < a < o"cr' ) and no extrema for a > leads us to expect longitudinal instabilities 
in the region u^\a) < u < u^\a) of u (cf. Fig. Id) in the first case, and no 
instabilities (cf. Fig. lb) in the second case. 5 To verify that the lowest eigenvalue Uq 
does indeed change sign at these values, we performed a numerical analysis whose results 
are shown in Fig. 5b and indeed reproduce the expected behavior. 6 Finally, for the 
angular fluctuations, the relevant Schrodinger potentials are given by the complicated 
expressions in Eqs. (A. 13) of appendix A. Starting from V e , we find that for a < 0.71 
it is positive while for a > 0.71 it develops a negative part which means that it can in 
principle support a bound state of negative energy. The critical value aiP above which 
5 For a > 1/4 the derivative term in (3.14) has definite sign, so (3.14) has no solution for any value 

of Mq. 

6 For this case, we can use a Schrodinger description for the longitudinal fluctuations and apply our 

(v) 

infinite- well approximation. This docs indeed lead to two critical values of u for < a < <j„ ■ 



a 

E~-L . 



(4.43) 
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Figure 5: (a) Stability diagram in the (a, u ) plane for the deformed sphere and the 
trajectory with (6o,ipo) — (0, 7r/4). The three shaded regions correspond to instabilities 
due to longitudinal, 66 and Sip fluctuations (lightest to darkest). The origin (0, 1) actually 
belongs to the stable part of the diagram. Specifically it is stable on its own right, 
but nowhere in its vicinity (for an analogous explicit example see sec. 5.4 of [13]). 
(b) Evolution of the lowest eigenvalue of the longitudinal fluctuations with u , plotted 
for a — 0, 0.15 < criP and 0.5 > a& (bottom to top), (c) Same for the angular 59 
fluctuations, for a = < a%P , 1 > o'er"* and 10 3> a%P (top to bottom). 

such states can exist is determined in the infinite-well approximation to be a^P ~ 0.78, 

( 0) 

which is quite close to the true value <7cr — 0.805 determined by numerical analysis. As 

(0) 

a increases beyond this value, the corresponding critical value u 0c approaches the value 
of 1.097 as a becomes very large. Turning to the Sip fluctuations, the same reasoning 
used in the deformation of the conformal background for the same trajectory carries over 
to the present case, implying that there exists a critical value affl = 1/ a/15 (within our 
numerical limitations) above which there occur instabilities for all values of uq. This value 
is presumably inherited by the conformal limit of the potential as analyzed in subsection 
4.3.1. The above results are summarized in the stability diagram of Fig. 5. 

We also note that, in the limit a >> 1, the potential for the 86 fluctuations reads 



and x can be explicitly determined in terms of u by (4.53). In this limit, the problem 
can be solved exactly (see appendix B), leading to the critical value u^) ~ 1.09737, by 
mapping it to the well studied Lame equation. 

From the above we conclude that for a finite narrow range of values of the deformation 
parameter, namely for 0.209 ~ < <j < criP — 0.258 the linear confining behavior 
for large distances is stable, after which it is wiped out first by the instability of the Sip 



V e = 2{u 2 - 2) + 0{a- 2 ) , 



(4.44) 
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fluctuations and subsequently by the 89 ones. Moreover, even for smaller values of a, 
namely for < a < — 0.209, the confining behavior of potential is stable as we have 
discussed already. 

• (0 O = 0,^0 = or 7r/2). The length and potential energy are the same as in the 
corresponding undeformed case, namely they are given by [4] 

L(u ) = 4^- [n(k'\ k) - K(k)} , (4.45) 

and 

E(uo) = ^2u 2 -l [k' 2 K(k) - E(k)] + 1 , (4.46) 

where 

k = ; U ° , k' = VT^¥ . (4.47) 
Their behavior is as in Fig. lc. 

By our general results, there exist longitudinal instabilities occurring for u below the 
critical value for the undeformed case, namely uq c ~ 1.125. For the angular 89 and 8ip 
fluctuations, the Schrodinger potentials are given by (A. 14) of appendix A and differ by 
the ones in the undeformed case by the positive-definite term 4cr 2 « 2 , while the expression 
for the Schrodinger variable x in terms of u is the same as in the undeformed case. The 
above results imply that, since the 89 and 8ip fluctuations are stable in the undeformed 
case, they are stable in the present case as well. 

• (# ,^o) = (7r/2,any). In this case, the length and potential energy are again the same 
as in the undeformed case and are given by [4] 



L ( u o) = — 



n(i k)-K(k) 



(4.48) 



and 



where 



E{uo) = ^ [K(k) - 2E(k)] . (4.49) 



k = J^l, k' = VT^W. (4.50) 
V 2u 

The behavior is as shown in Fig. lb. For uq ^> 1 (small L), the behavior is Coulombic, 
whereas in the opposite limit, u — > 1 (large L), the asymptotics of (4.48) and (4.49) lead 
to the linear potential 

E ~ ^ , for L > 1 . (4.51) 
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In the undeformed case, the appearance of a linear confining potential is quite unexpected 
from the gauge-theory side and the paradox is resolved by the stability analysis of [13] 
which indicates that the configurations where this behavior appears are unstable under 
angular perturbations. However, for the deformed case, where M = 4 supersymmetry is 
broken to J\f = 1, the appearance of a confining potential is expected. Remarkably, this 
is confirmed by the stability analysis that follows: 

We first note that the transverse and longitudinal fluctuations are manifestly stable, the 
latter fact following from the absence of extrema of L{uq). For the angular fluctuations, 
the Schrodinger potential reads 

V e = -l + 4a 2 (u 2 -l) , (4.52) 

whence we see that the effect of the deformation is to modify V e (which, for the unde- 
formed case equals —1) by a positive-semidefinite term. It is quite intuitive that the 
addition of this term will tend to stabilize the angular fluctuations. To see how this 
occurs, we note that in the present case, the Schrodinger variable x can be explicitly 
determined in terms of u as 

where F(z/, k) is the incomplete elliptic integral of the first kind and k is the modulus 
defined in (4.50). Likewise, the endpoint x is determined in terms of u by 

hk 2 — 1 

Xq = x{uq) = y K(/c) . (4.54) 

Then, as shown in detail in appendix B, the equation for the angular fluctuations takes 
the form (B.5) with A and h given by the second of (B.7). For general values of a, this 
equation can be solved only numerically and the results are as shown in Fig. 6. We verify 
that the critical point Woe below which instabilities occur is a decreasing function of a 
which implies that the region where confinement occurs is gradually stabilized as one 
increases a. 

We finally note that for the special values 4<r 2 = n(n + 1), where n is a positive integer, 
(B.5) reduces to the Lame equation and can be solved exactly, as done in appendix B 
for the particular case n — 1. The resulting value of Uq c , given in the second of (B.14), 
agrees with our numerical results. 
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(a) (b) 

Figure 6: (a) Stability diagram in the (cr, tt ) plane for the deformed sphere and the tra- 
jectory with (8q,iPo) = (7r/2, any). The shaded region corresponds to instabilities under 
angular perturbations, (b) Evolution of the lowest eigenvalue of the angular fluctuations 
with u , plotted for a = 0, 0.25 and 0.5 (bottom to top). 

4.3.3 The disc 



We finally consider the deformation of the background corresponding to D3-branes dis- 
tributed on a disc of radius r , which is obtained from the deformed sphere background 
by the analytic continuation r 2 . — > — r%. The possible trajectories are the same as for the 
sphere. We examine them in turn: 

• (#0,^0) = (0, 7r/4). The length and potential energy are given by (see sec. 6.2 of [31]) 

2u 



L(u ,a) 



[n(a 2 ,A;)-K(A;)] , 



(4.55) 



1+^ 



and 



E(u , a) = VTT^ I ^2^ + _i_ [ a * K (k) - E(k)] + E(c) - K(c) 



where 



k 



Ur 



ul + 



l+a 2 



+ 



; Va- V 2m + 1 f „i 



1 ' 



C = 



(4.56) 



(4.57) 



For u 3> 1 (small L) we recover the standard Coulombic behavior enhanced by the factor 
\/l + a 2 , while for uq — > the asymptotics lead to the potential 

(7T - L) 2 



8E(kr) 



(4.58) 
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which shows that there is complete screening at the screening length 

L c = 7i , (4.59) 

that is invariant under a-deformations. 

Regarding stability, the solution is stable under transverse and longitudinal perturba- 
tions. Turning to angular fluctuations, the relevant Schrodinger potentials are given by 
the complicated expressions in Eqs. (A. 16) of appendix A. Starting from the 59 fluc- 
tuations, we use the infinite-well approximation to find that the bottom of the well is 
raised with increasing a and that there is a critical value > 1/ vTo" above which no 
instabilities exist. The existence of this critical value can be demonstrated by noting that 
for a >> 1, the potential takes the form 

V e ~ 2{u 2 + 2) + 0(a' 2 ) . (4.60) 

This potential is positive and, using the techniques of appendix B, one can then show 
that perturbation theory holds about a — > oo. Thus, no instabilities occur in this limit. 
For the 5ip fluctuations, the same arguments used for the deformed conformal background 
for the same trajectory show that there exists another critical value atP = l/y/TE (again 
within numerical limitations and as before inherited by the conformal limit behavior of 
the solution) above which there occur instabilities for all values of uq. 

• (9 ,ip ) = (0,0 or 7r/2). In this case we find that the classical solution is the same as 
in the undeformed disc at 9 = [4] 

LW = M[n(^fc)-KW] (4.61) 

Uq -h 1 

and 

E(u ) = ^2u 2 + l [k' 2 K(k) - E(fc)] , (4.62) 

where 

k = , U ° , k' = , (4.63) 

Their behavior is as in Fig. la. 

By our general results, the solution is stable under transverse and longitudinal pertur- 
bations. For the angular 59 and Si/j fluctuations, the Schrodinger potentials are given 
by (A. 17) of appendix A and again differ by the ones in the undeformed case by 4a 2 u 2 , 
while the expression for the Schrodinger variable x in terms of u is the same. For the 69 
fluctuations, which have instabilities in the undeformed case, the extra term stabilizes 
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(0) 

the fluctuations and there exists critical value <7c r above which there are no instabilities. 
Finally, the Sip fluctuations are stable in the undeformed case and hence are stable in the 
present case as well. 

• (9 0: ip ) = (7i"/2,any). The length and potential energy are given by [4] 



2v 2 



ul 



n -jxr.fe - K ( fc ) 



(4.64) 



and 



where now 



E(u ) = t=^== K(k) - sjul + ul E(fc) + 1 , (4.65) 



'ul — u 2 



k=J > 7 , A;' = v / T^ (4.66) 
V M o + 

and (m<) denotes the larger (smaller) between w and 1. Their behavior is as in Fig. la 
and, in particular, we have a screened Coulomb potential with screening length 

L c = \ . (4.67) 

By familiar arguments, the solution is stable under transverse and longitudinal per- 
turbations, whose behavior is insensitive to the deformation. Regarding the angular 
fluctuations, the Schrodinger potential (equal to 1 in the undeformed case) now reads 

V e = l + 4a 2 (w 2 + l) , (4.68) 

and hence the solution is still stable under angular perturbations. 



5 Discussion and concluding remarks 

In this paper, we completed the analysis of our earlier work [13] by developing a formal- 
ism for studying the perturbative stability of strings dual to quark-antiquark pairs in a 
wide class of backgrounds with non-diagonal metrics and possibly a £>-field turned on. 
In particular, we derived a set of results by means of which the regions where instabilities 
occur can be determined by studying the zero modes of the differential equations govern- 
ing the fluctuations. The simplifying fact making this extension possible is that, although 
the various fluctuations are generally coupled, their zero modes actually decouple so that 
the problem can be treated in an analytic manner similar to the diagonal case. 
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The methods developed here were applied to strings in boosted, spinning and marginally- 
deformed D3-brane backgrounds, dual to quarkonium states. For the case of boosted 
D3-branes, we easily recovered the result that instabilities appear for the energetically 
unfavored "long" strings stretching further into the radial direction. For the case of 
spinning branes, we found regions of instabilities due to both longitudinal and angular 
fluctuations. Finally, for the case of marginally-deformed D3-branes, we found that 
angular fluctuations tend to completely destabilize the classical solutions for large enough 
values of the deformation parameter a, except for some special cases where the classical 
solutions are unaffected by the deformation and for which the deformation parameter a 
may actually have a stabilizing effect. In particular, we found parametric regions giving 
rise to a linear confining potential for which the dual string configurations are stable. 

As the methods developed here are completely general, we expect that they are di- 
rectly applicable to other backgrounds and, in particular to other gravity duals of gauge 
theories with reduced or no supersymmetry. As such backgrounds typically involve a 
non-trivial dependence on certain angular coordinates, it would be particularly interest- 
ing to determine the parametric regions for which quarkonium string configurations in 
these backgrounds are stable under angular perturbations. Extending our formalism to 
fluctuations of general brane probes should be also possible. 
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A Angular Schrodinger potentials 

Here we collect, for reference purposes, the Schrodinger potentials for angular fluctuations 
for all metrics and trajectories considered in the present paper as well as in [13]. 
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Multicenter D3-branes 

The only "angular" coordinate is 9. The Schrodinger potentials for the various trajecto- 
ries are given below. 

• Sphere, 9 = 0: 

2u* -u* + kk - 1) - iK - jjugQig - 1) 
Mo) 4^ _ 1)2 • ( A - X ) 

• Sphere, # = tt/2: 

^(«;«o) = -l. (A.2) 

• Disc, #o = 0: 

2 M 6 + u 4 + [GijK + 1) ~ IK + + 1) 

^(«; ««) = 4u2(U 2 + rp • 

• Disc, # = tt/2. 

V^(u;uo) = l. (A.4) 

Spinning D3-branes 

The only "angular" coordinate is 9 as before. The Schrodinger potentials are given below. 

• Two equal angular momenta, 9 = 0: 

V e (u;u ,X) = ] _ {2u w - u 8 + {6[ug(tig - 1) + A 4 ] - l}u° 

- [3u 2 (u 2 - 1) + 11AV - h[2ul(ul - 1) - 1]AV 
+ 7ul{ul-l)\*) . (A.5) 

• Two equal angular momenta, 9 = n/2: 

Ve(u;u ,X) = -1 . (A.6) 

• One angular momentum, 9 = 0: 

Ve(u;u ,X) = - * - {2u w + 3^ 8 + {6frK + 1) + A 4 ] - 1} M 6 

+ {9[«jj(ug + 1) + A 4 ] - 1}m 4 + [(3 - 10A 4 )ug(ug + 1) + 3A> 2 

- 5«g(«g + l)A 4 }. (A.7) 

• One angular momentum, 9 = n/2: 

V e (u;u ,X) = 1 • (A.8) 
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Deformed D3-branes 

Now, the "angular" coordinates are 9 and ip, with ip becoming irrelevant for the trajec- 
tories with 9 = 7r/2. The Schrodinger potentials are as follows. 

• Conformal, (9 ,ip ) = (0, 7r/4): 

2a 2 4<T 2 
V e (u; u , a) = — — -u 2 , V^u; u , a) = -— — -u 2 . (A.9) 
1 + a z 1 + a z 

• Conformal, (9 0: t/j ) = (0,0 or ir/2): 

Vo(u; u , a) = Vf(u; u , a) = Aa 2 u 2 . (A.10) 

• Conformal, (#0,^0) — ( 7r /2,any): 

Ve(u;u ,a) = Aa 2 u 2 . (A.ll) 

• Conformal, (9 ,ipo) = (sm'^l/VS), tt/4): 

8a 2 

Vo(u; u , a) = V^(u; u , a) = - - ^ 2 u 2 . (A.12) 

• Sphere, (0 O ,^) = (0,vr/4): 

V e (u;u ,a) = ^ - _ ^ {8a 2 (l + a 2 )u 10 + 2(1 + a 2 )(l - 15a 2 + 8a> 8 

- (3 - 41a 2 - 36a 4 )w 6 + 3{2(1 + a 2 )w 2 [(l + a 2 )u 2 - 1] - 7a 2 }w 4 

- {(9 + 8a 2 )[(l + a 2 )u 2 - 1] - 1}m 2 + 3w 2 [(l + a 2 )u 2 - 1]} , 

^(«;«o.-) - - 4 , 2 (, 2 -!)[(! + .> 2 -l]3 { 16ff2(1 + - V2 (A ' 13) 

- 2(1 + a 2 )(l + 35a 2 + 18<r 4 y° + (5 + 103a 2 + 94a 4 + 12a 6 )u 8 

+ {6(1 + a 2 )(l + 2a 2 )w 2 [(l + a 2 )u 2 - 1] - 3 - 54a 2 - 20a 4 }w 6 

- {(13 + 18a 2 + 4a 4 )« 2 [(l + a 2 )u 2 - 1] + 1 - 7a 2 }w 4 
+ { 8M 2 [(l+a 2 ) M 2 -l] + lK- M 2 [(l + a 2 K-l]} . 

• Sphere, (0 o ,^o) = (0, or tt/2): 

Vg(u; u , a) = V g (u; u ) + Aa 2 u 2 , V^{u\ u Q , a) = V^(u; u ) + 4aV , (A.14) 

where Vo(u;uq) and V^(u;uo) are the potentials for the corresponding undeformed case, 
with the former given by (A.l) and the latter given by a function which, since the ip 
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fluctuations in the multicenter case are stable, admits no negative-energy bound states. 

• Sphere, (9 ,ip ) = (7r/2,any): 

V g (u;u ,a) = -1 + 4a 2 (w 2 - 1) . (A.15) 

• Disc, (0 o ,Vo) = (0,tt/4): 

V e (u;u ,a) = ^ - 1^ - i]3 { 8 a 2 (l + <r>" - 2(1 + a 2 )(l - 15a 2 - 8aV 

- (3 - 41a 2 - 36a 4 )w 6 - 3{2(1 + a 2 )w 2 [(l + a 2 )w 2 + 1] - 7a 2 }« 4 

- {(9 + 8a 2 )[(l + a 2 )u 2 Q + 1] - 1}« 2 - 3« 2 [(1 + a 2 )w 2 - 1]} , 

= - i^T i)[(i + a vTf { 16ff2(1 + g>12 (A - 16) 

+ 2(1 + a 2 ) (1 + 35a 2 + 18<7 4 )u 10 + (5 + 103a 2 + 94a 4 + 12a 6 )w 8 

- {6(1 + a 2 )(l + 2a 2 )« 2 [(l + a 2 )u 2 + 1] - 3 - 54a 2 - 20a 4 }w 6 

- {(13 + 18a 2 + 4a 4 )w 2 [(l + a 2 )u 2 + 1] + 1 - 7a 2 }w 4 

- {8m 2 [(1 + a 2 )ul + 1] + 1}m 2 - m 2 [(1 + a 2 )ul + 1]} . 

• Disc, (0 o ,^o) = (0,0 or tt/2): 

V d (u; u , a) = Vg(u; u ) + Aa 2 u 2 , V^(u; u , a) = V^(u; u ) + 4aV , (A.17) 

where V$(u; u ) and V^(u; uq) are the potentials for the undeformed case, with the former 
given by (A. 3) and the latter by a function admitting no negative-energy bound states. 

• Disc, (0 O , ^0) = (tt/2, any): 

V e (u;u ,a) = l+4a 2 (« 2 + l) . (A.18) 



B Angular fluctuations and the Lame equation 

For the special cases of (i) the deformed sphere at O = in the limit a 3> 1 and (ii) 
the deformed sphere at 0o = 7r/2 (for any a), the Schrodinger equation for the angular 
fluctuations takes a relatively simple form that allows for explicit solutions. In both 
cases, the Schrodinger variable x is given in closed form in terms of u by Eq. (4.53), 
which can be inverted for u to yield 

" 2 = 2 < Uvw) " 9 = 2< { p( ^ x) + is) • 
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where sn(z, k) is the elliptic Jacobi function and p(z) = p(z \ 001,002) = p(z,g 2 ,g 3 ) is the 
Weierstrass p-function, specified by the half-periods 



1 +v 2 

wi = K(fc), ou 2 = iK(k>), fc* = -ZJo, (B.2) 



or, equivalently, by the elliptic invariants 

g 2 = -4(eie 2 + e 2 e 3 + e 3 ei) , g 3 = 4eie 2 e 3 , (B.3) 

with the ej given by 

3u§-l 1 3^ + 1 

ei = " 2 , e 2 = —J , e 3 = , B.4 

6-Uq 3-Uq 6-Uq 

and satisfying ei + e 2 + e 3 = 0. Substituting (B.l) into (4.44) and (4.52) to obtain explicit 
expressions for the potentials in terms of x, and making a further change of variables to 
z = \^2uqx, we arrive at the following equation 

d 2 1 oj 2 

+ [Ap(z) + h]}V(z) = —zV(z), 0O<wi, (B.5) 



dz 2 ' L " rK ~' ' " J J 2u\ 
subject to the boundary conditions 

tf(0) = , = , (B.6) 

and to the requirement that ^(2) be real. Here, A and ft, stand for the constants 

5 



(i) : A = 2 , h 



3u 2 ' 



(ii) : A = 4a 2 , h = ■ (B.7) 

om 



In the special cases where A = n(n + 1) (with n being a positive integer) Eq. (B.5) is 
just the Lame equation, whose exact solutions are well-known (see e.g. [38]). Although 
the construction of the specific solutions satisfying the boundary conditions (B.6) is 
rather cumbersome, it is instructive to consider the simplest possible case, n — 1, as a 
consistency check of our numerics. For this case, the zero-mode equation simplifies to 

d 2 1 

+ [2p(z) + h] \V(z) = , O^z^Ui. (B.8) 



dz 2 

For the analysis that follows, it is convenient to define the quantities 

a = p^ih) , oo 3 = ui + oj 2 . (B.9) 
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When a j£ u>i, i = 1,2,3, the equation (B.8) possesses the linearly-independent solutions 

* ± ( 2! )=exp[^C(a)]^^ , (B.IO) 

where C(z) = C(z \ cu 1 , u 2 ) and u(z) = u(z \ u>i, o; 2 ) are the Weierstrass zeta and sigma 
functions. 7 When a = Ui, the linearly-independent solutions are instead 

* itl (z) = exp [-ztiwi)] ^i±^l , y i2 ( z ) = [Z(z + ui) + ze,} ^(z) . (B.ll) 

a{z) 

For these two cases, the general solution reads ^f(z) = C + ^ + (z) + C-^-(z) and ^(z) = 
Ci,i^i,i(z) + Ci^i^iz) respectively. Imposing the condition \I/(0) = 0, we obtain C + = 
C- and Ci t \ = — £(0^)6^2 respectively, while using the reality condition for ^(z) we find 
that a must be real. For the case a ^ Ui, the condition ^'(cui) = leads to the equation 

{a(cu! - a) [C(a) - CM + - a)]} {exp [2 (aCM - wiC(a))] - 1} = , (B.12) 

which is to be solved for a in the square {0 ^ Rea < 2uj 1 , ^ Ima < 2|w 2 |} m the 
complex plane, due to the periodicity of the Weierstrass functions under a — > a + 2mu 1 + 
2nu 2 - The trivial zeros of (B.12) occur at the points a — oj^ % — 1, 2, 3 where either 
the first or the second factor in curly brackets vanishes. However, for these points, we 
actually have to use the second set of linearly-independent solutions (B.ll). Doing so, 
we find that the values a = cj 2 ,3 are not allowed due to the reality condition, while for 
the value a = uj\ the boundary condition ^'{ui) = leads to the equation 3C(^i) = eii^i, 
which has no solution. Nontrivial roots of Eq. (B.12) can only arise from the second 
factor in curly brackets, i.e. they are given by the solutions of the equation 

aCM = wiC(a) , (B.13) 

subject to the condition a ^ uoi. Using Eqs. (B.7), (B.2) and (B.9), we can express this 
as a transcendental equation for u 0c , whose solutions for the two cases of interest are 

(i) : u 0c ~ 1.09737 , L c = 1.9312 , 

(ii) : M 0c ^ 1.02676 , L c = 2.92397 , for 4<r 2 = 2 . (B.14) 

Both results are in precise agreement with our numerical analysis. 
7 Some properties of these functions, used in the derivation of (B.12), are 

C'(z) = -p(z) , &{z) = o{z)t{z) , 

C(z + 2wi, 2 ) = C,(z) + 2C(wi, 2 ) , (t(z + 2wi, 2 ) = - exp[2(z + wi >2 )C(wi,2)M*0 • 
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C Soap films and strings in Rindler space 



An example of the calculations considered here and in [13] has been presented recently in 
[39] and refers to open strings ending on branes in Rindler space and their stability prop- 
erties. As is well-known, Rindler space is (a portion of) flat Minkowski space, expressed 
in coordinates adapted to an observer at constant acceleration. Such an observer experi- 
ences the Unruh effect, perceiving the inertial vacuum as a state populated by a thermal 
distribution of particles at a temperature T = k/2tt, where k is the surface gravity. It is 
then natural to expect that an open string with fixed endpoints propagating in Rindler 
space will have similar properties as one propagating in a black-hole background and, in 
particular, will exhibit a phase structure of the type shown in Fig. lc. In this appendix we 
prove the amusing fact that this problem is exactly equivalent to the classical-mechanical 
problem of the shape of a soap film suspended between two circular rings. 

The metric for Rindler space has the form 

ds 2 = -n 2 u 2 dt 2 + dy 2 + du 2 + ... , (C.l) 

where u is the radial direction with the Rindler horizon corresponding to u — and 
y is a generic spatial direction. The string configuration of interest corresponds to an 
open string with its two endpoints located at the same radius u = A and separated by 
a distance L along the y direction. It is obvious that the formalism of section 2 readily 
applies to this setup. Passing to convenient dimensionless units through the rescalings 

kA 2 

u — > Au , u -> Au , L — > AL , E — > — -E , (C.2) 

An 2 

inserting the metric components of (C.l) in (2.12), and integrating from Uq to A, we find 
that the classical solution is the catenary curve of Leibniz, Huygens and Bernoulli, 

w = -u cosh— , (C.3) 

which is a slice of Euler's catenoid. The integration constant u and the energy of the 
string are determined by Eqs. (2.14) and (2.15) (without the subtraction term) which 
for our case read 

L = 2u [ . dU = 2M cosh~ 1 — , (C.4) 
a/m 2 -ul u 



■ii a 



and 

E = 4n I' ^£= = 2n( x fT^u 2 + u 2 cosh- 1 -) . (C.5) 



UQ 



- u 2 V V % 
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respectively. Eqs. (C.3)-(C.5) are exactly the same formulas appearing in Plateau's 
problem for a thin soap film of mass density 1 and surface tension | stretched between 
two coaxial rings of radius 1 that are separated by a distance L, with uq being the minimal 
radius reached by the film (see Eqs. (A.4)-(A.5) in appendix A of [13]), and thus the 
properties of the solutions (see e.g. [40]) are the same. Namely, there exists a critical 
value uo c of u and a corresponding value L c of L at which L'(u ) = 0, i.e. 

1 



1 -M^cosh" 1 — = 1 . (C.6) 

Solving this equation, we find Uq c ~ 0.552, whence L c ~ 1.325 (to compare with the 
results of [39], set u — > L/h). For separations L < L c , there is a "short" string (dual to 
the shallow catenoid), a "long" string of higher energy (dual to the deep catenoid), and a 
configuration of two straight strings (dual to the Goldschmidt solution). For separations 
L > L c , only the last configuration is possible. 

Regarding stability, the only relevant fluctuations are the longitudinal fluctuations along 
y. By the results of section 2, instabilities occur for u below the critical value where 
L'(uq) = i.e. for u < uq c . This can be verified directly by noting that the longitudinal 
Sturm-Liouville equation has the form 

(u 2 - ul)5y\u) + ^^5y\u) = -- 2 5y{u) . (C.7) 

The zero-mode solution satisfying the boundary condition 5y(l) = is given by 

5y{u) ~ cosh" 1 — M - cosh" 1 — + 1 , (C.8) 

u ^u 2 - ul Uo yjl - u 2 q 

and imposing the boundary condition lim u _^ u + [5y(u) + 2{u — u )Sy'(u)] = leads indeed 
to Eq. (C.6). Alternatively, one may set x = cosh -1 ^ and ^ = tanhx Sy to transform 
(C.7) to a Schrodinger equation with the Poschl-Teller type II potential V = 

COSil X 

and the energy defined in the interval x E [l,cosh _1 ^] and obeying Neumann and 
Dirichlet boundary conditions at the left and right endpoint respectively. 8 The solutions 
are the associated Legendre functions P^(tanhx) and <5i(tanhx), where v — — , and 
imposing the boundary conditions on the zero-mode solutions leads again to Eq. (C.6). 
In both formulations, it is very straightforward to check that the lowest eigenvalue really 
changes sign from positive to negative as uq crosses uq c from the right. In conclusion, the 
long strings with uq < u Qc are perturbatively unstable in the region where they exist. 



8 This is different than the corresponding Schrodinger equation for the fluctuations of the thin soap 
film problem (sec [13, 42]). 
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The qualitative analogy of the soap- film problem with AdS / CFT calculations of Wilson 
loops has already been considered in previous work [13, 41]. It is a pleasant surprise that 
this problem actually has an exact analog in the context of the Nambu-Goto string. 
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